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We use QCD sum rule approach to calculate the masses of the ground-state Aq and E^' baryons. 
Contributions of the operators up to dimension six are included in operator product expansion. 
The resulting heavy baryonic masses from the calculations are m^ h — 5.69 ± 0.13 GeV, and m^ a = 
2.31±0.19 GeV for A Q ; m Eb = 5.73±0.21 GeV, m s . = 5.81±0.19 GeV, m Ec = 2.40±0.31 GeV and 
m s j = 2.56 ± 0.24 GeV for Eq , respectively, which are in good agreement with the experimental 
qq ' values. 
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Recently, the CDF Collaboration has reported the observations of the heavy baryons Eb and [l| . Up 
to now, the masses of the ground-state Aq and Eq/| (Q = b, c) baryons as well as masses of several lower 
lying excited heavy baryons have been measured y, Q] . With the accumulation of experimental data, 



in 

O i 1 reliable and comprehensive theoretical explanations are needed. There were some theoretical investigations 
q^, on the heavy baryon masses, such as quark models [5, 6], mass formulas Q, and lattice QCD calculations 



[8|. With QCD sum rules |9j, heavy baryons were first discussed in heavy-quark limit in Ref. [10|, then 



masses for heavy baryons were calculated in the heavy quark effective theory (HQET) to leading and 
next-to- leading order in a s 11|, and to order 1/rriQ |l2j, Ilia]. In Refs. 14 1. the calculations for the heavy 



baryons began with the full theory and results of the calculations were exp anded by heavy-quark masses. 
Lately, the masses of Sq and H.q have been studied in QCD sum rules 15 , 3| . Along with the significant 



observation of Eb and Ej* , several renewed theoretical studies have been done by using various approaches 
mentioned above [l3| ■ In this paper we shall study heavy baryonic two-point correlators and obtain mass 
sum rules for Aq and Eq' 1 , using the technique developed in 3 IS I- The paper is organized as follows. 



In Sec |TT] we have derived QCD sum rules for Aq and Eq . Section iHll is devoted to numerical analysis 



and discussions. This section also includes a brief summary. 



(*) 



II. QCD SUM RULES FOR Aq AND E£ 

The basic points in the application of QCD sum rules to problems involving heavy baryons are to choose 
certain suitable interpolating currents in terms of quark fields. The currents for Aq and doublet {Eq, Eq} 
are associated with the spin-parity quantum numbers j p = + and j p = 1 + for the light diquark system 
with antisymmetric and symmetric flavor structure, respectively. Adding the heavy quark to the light- 
quark system, one obtains j' p = l/2 + for the baryons Aq and the pair of degenerate states j p = 1/2+ 
and j p = 3/2 + for the baryons £q and Eq. Consequently, we need an isospin-0 quark pair for Aq and 
an isospin-1 quark pair for Eq, after adding a third quark to form antisymmetric and symmetric flavor 
structure respectively, which can determine the choice of 7 matrices in baryonic currents [Io| . For Eg, 



the currents can be obtained from those of Eq using SU(3) symmetry relations 20(. Concretely, we adopt 
le followinj 
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the following forms of currents for the heavy baryons Aq and Eg , which can be generally written as 



JQ = s abc (qI a cr k q b 2 )r' k Q c 



(1) 
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for jP = 1/2+, and 

j Q = e abc [2/V3(qf a Cr k Q b )f k q^ + 1/V%f a CT k q\)T k Q c ] (2) 
for j p — 3/2+ . Tfc and T fc are chosen covariantly as 

r fc = 75,r; = i (3) 

for Aq baryons, and 

r fc = 7^ r fc = 7^75 (4) 

for Eg baryons. Here the index T means matrix transposition, C is the charge conjugation matrix, and a, 
b, c are color indices. The QCD sum rules for Aq and Eg are constructed from the two-point correlation 
function 

n(g 2 ) = i J dW^QlTfoCaOJgCOjp). (5) 

Lorentz covariance, parity, and time reversal imply that the two-point correlation function in Eq. ((3J) has 
the form 

U(q 2 )=Il 1 (q 2 )+4IL 2 ( q 2 ). (6) 

According to the philosophy of QCD sum rules, the correlator is evaluated in two ways. Phenomenologi- 
cally, the correlator can be expressed as a dispersion integral over a physical spectral function 

, 2 , . 2 i + M H 1 f°° , ImnP hen (s) 

n(g 2 = X% f-5 + - / ds + subtractions, 7 

Mfj -q 2 7T J so s- q 2 

where Mjj denotes the mass of the heavy baryon. In the operator product expansion (OPE) side, short- 
distance effects are taken care of by Wilson coefficients, while long-distance confinement effects are included 
as power corrections and parameterized in terms of vacuum expectation values of local operators, the so- 
called condensates. Hence 

n^ 2 ) = nfV) + nf ond (g 2 ), i = 1,2. (8) 

We work at leading order in a s and consider condensates up to dimension six. To keep the heavy-quark 
mass finite, we use the momentum-space expression for the heavy-quark propagator. We follow Refs. 
and calculate the light-quark part of the correlation function in the coordinate space, which is then 



Fourier-transformed to the momentum space in D dimension. The resulting light-quark part is combined 
with the heavy-quark part before it is dimensionally regularized at D = 4. For the heavy-quark propagator 
with two and three gluons attached we use the momentum-space expressions given in Ref. 22| . With Eq. 
© , we can write the correlation function in the OPE side in terms of a dispersion relation 



n i (g a )=/ d s -^ + n™ d ( q % (9) 

s q 



where the spectral density is given by the imaginary part of the correlation function 

Pi (s) - -Imnp PE ( S ). (10) 

Equating the two expressions for Tl(q 2 ) and assuming quark-hadron duality yield the sum rules, from 
which masses of the heavy baryons can be determined. After making a Borel transform and transferring 
the continuum contribution to the OPE side, the sum rules can be written as 

\\M H e- M »l M * = r ds Pl (s)e- s / M2 + Bnf ond (11) 

Jm 2 Q 
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\ 2 H e- M *' M2 = / d S p 2 (s)e-^ M2 + M™ nd . (12) 

To eliminate the baryon coupling constant Xh and extract the resonance mass Mh, we first take the 
derivative of Eq. (|11[) with respect to 1/M 2 , divide the result by itself and deal with Eq. (fT2"|) in the same 
way to get 

Ml = { P d Spl {s)se-»/ M2 + d/d(- -L )5IT™ nd ( s )}/{ f° d S/5l ( S )e-^ 2 + MI? ond ( s )} (13) 
M| = {/ ^ 2 ( S ) S e- s / M2 +d/rf(--l)Mr d W}/{ d S p 2 ( s ) e -/M 2 + snf ond ( s )}, (14) 

./nA M J m 2 



where 



with 



for Aq baryons, 



Pi ( S )= P r rt (s)+p^(s)+pf\ S ) as) 



-^mg I da{ -f{m\ \- saf (16) 

^00 = ^- Q / 1 ^[(^) 2 + 2] (17) 



P? ert W = f daV-^{ml-saf (19) 

/r 2 \ (q 2 G 2 ) ml, „ 

4 G> « = ^Mi- ^ 2 (20) 

M cond = _i9^) m 2 /' 1 da( 1 _^ ) 2 e -^/( Q M 2 ) 
2 3 • 2 9 tt 4 Q J Q a 7 



(21) 



P? ert (s) = 2^I^Q J\ a ^f{ml-saf (22) 
/^M = (23) 
3 • 2 t> 7r 4 * Jo ar 
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Pr\ S ) = ^fda^L^-saf (25) 

Z 7T ./ a Cx 



A 

pf\s) = -^[l- {m y s f] (26) 



2 7 7T 



/ill?"" 1 = - j^™ 2 Q f ' da{ — fe <• 
, 4(gg) 2 



o 



3 



a 

e -m 2 Q /M' 2 



da (t^ifn _ ^ )e -m 2 «/(^) (27) 

a aM z 



for Eg baryons, and 



^ 6rt(s) = 2^™° / ^(^) 2 K- S a) 2 (28) 



A 



a 



Pi q) ( s ) = 3^(59) f daiml-^ ) (29) 



A 



= ^ W( =/ A ^ (1 ^ )2 - 2] (30) 



o n COnd _ (ffi 2 ) 3 / , (I-")' VJIaM') 



pert,, 1 (l-a) 2 (-l + 3a) , 2 2 
^2 0) = orZJ / da ^2 ( m Q- sa ) ( 32 ) 



2 5 7T 

P^OO = -^(qq)m Q (l-m%/s) 2 (33) 

{G 2 } / \ (9 2 G 2 ) 1 , (my S ) 2 /- 1 2 
P 2 ( s ) = 3.26^4 [~2 + — 2 / ' ~ >( m Q ~ sa >* ( 3 > 

Wl2 - 3 2 -2v m «y Q rfa 

3tt z 

■ 4 <W> 2 e - m 2 Q /M 







3 2 

<-9 3G3 ) jf* dai^[ a (l - 4a - 3a 2 ) - 2(1 - 4a + a 2 )^f]e-™^ M2 > (35) 



3 2 -2 9 7r 4 y a 
for Eg baryons. The lower limit of integration is given by A = rrVq/s. 

III. NUMERICAL RESULTS AND DISCUSSIONS 

In the numerical analysis of the sum rules obtained above, the input values used for the quark masses and 
condensates are taken as: m c = 1.25 ± 0.09 GeV, m b = 4.20 ± 0.07 GeV [23] with (qq) = -(0.23) 3 GeV 3 , 
(gqa-Gq) = m 2 (gg),m 2 = 0.8 GeV 2 , (g 2 G 2 ) = 0.5 GeV 4 , and (g 3 G 3 ) = 0.045 GeV 6 jig. The proper areas 
of the thresholds can be determined from the consideration that the stability of the Borel curves should not 
be sensitive to them. According to the standard criterion in QCD sum rules, the Borel windows are fixed 
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FIG. 1: In (a) the dashed line shows the relative pole contribution (the pole contribution divided by the total, 
pole plus continuum contribution) and the solid line shows the relative continuum contribution from Eq. (12) for 
y^io = 6.5 GeV for £5. Its OPE convergence is shown in (b) by comparing the perturbative, quark condensate, 
two-gluon condensate and three-gluon condensate contributions. 

in such a way 0,0,13 

: on one hand, in comparison with the condensate contributions, the perturbative 
contribution should be larger, and the lower limit constraint for M 2 in the sum rule windows is obtained; on 
the other hand, the upper limit constraint is obtained by imposing the restriction that the QCD continuum 
contribution should be smaller than pole contribution. Giving an illustration, the comparison between pole 
and continuum contributions from Eq. (12) for t/sq — 6.5 GeV for Ef, is shown in Fig. 1(a), and its OPE 
convergence by comparing the different contributions is shown in Fig. 1 (b) . The analysis for the others has 
similarly been done, but the corresponding figures are not listed in the paper for conciseness. Accordingly, 
the thresholds and Borel windows are taken as y/Jo = 6.2 — 6.4 GeV, M 2 — 4.5 — 6.0 GeV 2 for A&, 
= 2.7 - 2.9 GeV, M 2 = 1.5 - 3.0 GeV 2 for A c , y/sE = 6.4 - 6.6 GeV, M 2 = 4.5 - 6.0 GeV 2 for E b and 
EJ, = 3.0 - 3.2 GeV, M 2 = 1.5 - 3.0 GeV 2 for E c , and 3.1 - 3.3 GeV, M 2 = 1.5 - 3.0 GeV 2 

for E*. The Borel curves for the dependence on M 2 of the heavy baryon masses are shown in Figs. 2-4. 

In Table HI we present our results for the masses of Aq and Eg baryons and a comparison with 
experimental data and other theoretical approaches. In order to decrease the systematic errors of the sum 
rules we take the average of the results obtained from sum rules (|1 3[) and (|14p in the numerical evaluation. 
The errors reflect the uncertainty due to sum rule windows only; the uncertainty due to the variation of 
the quark masses and QCD parameters is not included. It is worth noting that the QCD 0(a s ) corrections 
in the perturbative expansion of the OPE have not been included in the sum rule calculations. However, 
it is expected that the QCD 0(a s ) corrections might be under control since a partial cancellation occurs 
in the ratio obtaining the mass sum rules (fTB")) and (fT4")) . This has been proved to be true in the analysis 



for the heavy baryons in the HQET [llj and for the heavy mesons in full QCD 22|, l25j . 



In summary, we have applied the QCD sum rule approach to calculate the masses of the heavy baryons 
Aq and Eg including the contributions of the operators up to dimension six in OPE. The final results 
extracted from our sum rules are: m^ b — 5.69±0.13 GeV, m\ c = 2.31±0.19 GeV, ms b = 5.73±0.21 GeV, 
m s * = 5.81 ± 0.19 GeV, m Sc = 2.40 ± 0.31 GeV, and m s . = 2.56 ± 0.24 GeV. The gained masses are 
well compatible with recent experimental data. However, there are still some differences from our central 
values and experimental values, which implies that the predictions might be improved by the computation 
of the QCD 0(a s ) corrections. 
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FIG. 2: The dependence on M for the masses of A;, and A c . The continuum thresholds are taken as y/so = 
6.2 — 6.4 GeV, y^So = 2.7 — 2.9 GeV respectively: (a) and (b) are from sum rule (|13p . (c) and (d) from sum rule 
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FIG. 3: The dependence on M 2 for the masses of £&, EJ. The continuum thresholds are taken as ^/~sq = 6.4 — 
6.6 GeV: (a) and (b) are from sum rule (I13[l . (c) and (d) from sum rule (|14|l . 
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FIG. 4: The dependence on M 2 for the masses of E c and E^. The continuum thresholds are taken as y^so = 
3.0 — 3.2 GeV, y^so = 3.1 — 3.3 GeV respectively: (a) and (b) are from sum rule (|13p . (c) and (d) from sum rule 



TABLE I: the mass spectra of Aq and Eq 
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Baryon 


Experiment (MeV) 




Theory 














This work (GeV) Ref. [5] (MeV) Ref. 


[6] (MeV) Ref. [7] (MeV) Ref. [8] (MeV) Ref. [14] (GeV) 


A„ 


5619.7 ± 1.2 [4] 
5624 ± 9 [23] 


5.69 ±0.13 


5622 


5585 


5620 


5672 






A c 


2286.46 ±0.14 [3, 23] 


2.31 ±0.19 


2297 


2265 


2285 


2290 






E 6 


5807.8^2:2 ± 1-7 for E+ [1] 
5815.2 ± 1.0 ± 1.7 for E" [1] 


5.73 ±0.21 


5805 


5795 


5820 


5847 


5.70 n 


- 6.62 


XI 


5829.0±i;|±iJ for EJ+ [1] 
5836.4 ± 2.0li;? for E t~ [U 


5.81 ±0.19 


5834 


5805 


5850 


5871 


5.4 n 


- 6.2 


E c 


2453.76 ±0.18 [23J 


2.40 ±0.31 


2439 


2440 


2453 


2452 


2.45 <- 


- 2.94 


E c 


2518.0 ±0.5 [23] 


2.56 ±0.24 


2518 


2495 


2520 


2538 


2.15 - 


- 2.92 
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